#include "tbparameters.hpp"

struct tb parameters(string name,string mat){
  struct tb p;

  p.name  = name;
  p.mat   = mat;

  //-----------------------------------------------------------------
  if (name=="vogl"){
    if (mat=="gaas"){
      p.Esa    =-8.3431;
      p.Epa    = 1.0414;
      
      p.Esc    =-2.6569;
      p.Epc    = 3.6686;
      
      p.Essa   = 8.5914;
      p.Essc   = 6.7386;
      
      
      p.Vss    =-6.4513;
      p.Vxx    = 1.9546;
      p.Vxy    = 5.0779;
      p.Vsapc  = 4.4800;
      p.Vscpa  = 5.7839;
      p.Vpcsa  = p.Vsapc;
      p.Vpasc  = p.Vscpa;
      
      p.Vssapc = 4.8422;
      p.Vpassc = 4.8077;
      p.Vpcssa = p.Vssapc;
      p.Vsscpa = p.Vpassc;
      
      p.Vssss  = 0;
    }
    if (mat=="alas"){
      p.Esa    =-7.5273;
      p.Epa    = 0.9833;
      
      p.Esc    =-1.1627;
      p.Epc    = 3.5867;
      
      p.Essa   = 7.4033;
      p.Essc   = 6.7267;
      
      p.Vss    =-6.6642;
      p.Vxx    = 1.8780;
      p.Vxy    = 4.2919;
      p.Vsapc  = 5.1106;
      p.Vscpa  = 5.4965;
      p.Vpcsa  = p.Vsapc;
      p.Vpasc  = p.Vscpa;
      
      p.Vssapc = 4.5216;
      p.Vpassc = 4.9950;
      p.Vpcssa = p.Vssapc;
      p.Vsscpa = p.Vpassc;
      
      p.Vssss  = 0;
    }
    if (mat=="ge"){
      p.Esa    = -5.88;
      p.Epa    =  1.61;
      p.Essa   =  6.39;
      p.Vss    = -6.78;
      p.Vxx    =  1.61;
      p.Vxy    =  4.90;
      
      p.Vscpa  =  5.4649;
      p.Vssapc =  5.2191;
      p.Vpcsa  = p.Vsapc;
      p.Vpasc  = p.Vscpa;
      
      p.Esc    = p.Esa;
      p.Epc    = p.Epa;
      p.Essc   = p.Essa;
      p.Vsapc  = p.Vscpa;
      p.Vpassc = p.Vssapc;
      
      p.Vsscpa = p.Vpassc;
      p.Vpcssa = p.Vssapc;
      
      p.Vssss  = 0;
    }
  }
  if (name=="klimeck"){
    if (mat=="gaas"){
      p.Esa    = -3.53284;
      p.Epa    = 0.27772;
      p.Esc    = -8.11499;
      p.Epc    = 4.57341;
      p.Essa   = 12.33930;
      p.Essc   = 4.31241;
        
      p.Vss    = -6.87653;
      p.Vxx    =  1.33572;
      p.Vxy    = 5.07596;
        
      p.Vsapc  = 2.85929;
      p.Vscpa  = 11.09774;
      p.Vpcsa  = p.Vsapc;
      p.Vpasc  = p.Vscpa;
            
      p.Vssss  = 0;
      p.Vssapc = 6.31619;
      p.Vsscpa = 5.02335;

      p.Vpcssa  = p.Vssapc;
      p.Vpassc = p.Vsscpa;
        
      p.soa    = 0.32703/3;
      p.soc    = 0.12000/3;
      }
    if (mat=="alas"){
      p.Esa    =-3.1537;
      p.Epa    =-0.09711;
    
      p.Esc    =-9.52462;
      p.Epc    = 4.97139;
        
      p.Essa   = 12.05550;
      p.Essc   = 3.99445;
        
        
      p.Vss    =-8.84261;
      p.Vxx    =-0.01434;
      p.Vxy    = 4.25949;
      p.Vsapc  = 2.42476;
      p.Vscpa  =13.20317;
      p.Vpcsa  = p.Vsapc;
      p.Vpasc  = p.Vscpa;

      p.Vssapc = 5.83246;
      p.Vpassc = 4.60075;
      p.Vpcssa = p.Vssapc;
      p.Vsscpa = p.Vpassc;
    
      p.soa    = 0.29145/3;
      p.soc    = 0.03152/3;
    }
  }

  if (name=="chadi")
    if (mat=="ge"){
      p.Esa    = -5.85;
      p.Epa    =  1.50;
      p.Vss    = -6.75;
      p.Vxx    =  1.60;
      p.Vxy    =  5.40;
      p.Vscpa  =  5.33;
      p.soa    = 0.29/3;
      
      p.Esc    = p.Esa;
      p.Epc    = p.Epa;
      p.Essc   = p.Essa;
      p.Vsapc  = p.Vscpa;
      p.soc    = p.soa;
      
      p.Vpcsa  = p.Vsapc;
      p.Vpasc  = p.Vscpa;
    }

  return p;

};
